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Abstract. We study the implications of translation invariance on the tangent 
dynamics of extended dynamical systems, within a random matrix approxima- 
tion. In a model system, we show the existence of hydrodynamic modes in the 
slowly growing part of the Lyapunov spectrum, which are analogous to the hydro- 
dynamic modes discovered numerically by [Dellago, Ch., Posch, H.A., Hoover, 
W.G., Phys. Rev. E53, 1485 (1996)]. The hydrodynamic Lyapunov vectors loose 
the typical random structure and exhibit instead the structure of weakly perturbed 
coherent long wavelength waves. We show further that the amplitude of the per- 
turbations vanishes in the thermodynamic limit, and that the associated Lyapunov 
exponents are universal. 



1 Introduction 

The Lyapunov spectrum is recognized as an important diagnostic of chaotic dy- 
namical systems. As such, it has been studied intensely in the context of extended 
systems flR 0j. It has been observed that in the thermodynamic limit the spec- 
trum seems to approach a continuous density, and some theoretical studies have 
focused on this phenomenon |TJ]. However, despite the large amount of available 
data, there is an unsatisfactory degree of understanding of the relation between the 
Lyapunov spectrum of extended systems, and their global or collective properties. 

In connection with these problems, the recent study of [Q] presents an inter- 
esting development (see also [|5p for further results and references). In the con- 
text of molecular dynamics simulations, they find hydrodynamic, i.e., slow, long- 
wavelength behavior in the tangent space dynamics. Namely, they observe that 
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the Lyapunov vectors associated with the Lyapunov exponents of small absolute 
value have ordered, wave-like structure, and that the exponents themselves follow 
an ordered pattern. Hydrodynamic behavior in phase space is of course present 
in every extended system with a continuous symmetry. In the models considered 
in 0] the symmetries in question are translation and Galilei invariance, precisely 
those which enable the hydrodynamic description of a fluid in terms of the Navier- 
Stokes equations [j^] . However, it is for the first time that a similar phenomenon 
is observed in tangent space. 

In this paper, we study theoretically the slow Lyapunov modes (vectors and 
exponents) of extended systems with translation invariance. We focus attention 
on a simplified model which shares the essential features with the more elaborate 
model of ffi\. This simplified model is constructed only in tangent space without 
an accompanying real space dynamics, and is based on a random matrix approx- 
imation. As has been often found before, in systems with strong chaos, quali- 
tative features of the Lyapunov spectrum are well reproduced by approximating 
the tangent matrices by independent random matrices with appropriately chosen 
distributions [0, [5J]. We prove several statements on the slow Lyapunov modes of 
this model in the thermodynamic limit, which show that in this limit the Lyapunov 
vectors and exponents are indeed well described as being hydrodynamic. 

The basic reason for the existence of these hydrodynamic modes is evidently 
the translation invariance. Its presence dictates that the dynamics are indifferent 
to a uniform shift of all the particles (or their momenta), so that the associated 
Lyapunov vectors are decoupled from the rest of the dynamics, and the associ- 
ated Lyapunov exponents vanish. We show that slowly growing large wavelength 
disturbances are nearly decoupled for the same reason, and use this property to 
show how the clean wave structure is obtained as a result of the orthogonalization 
procedure which involves all the faster growing Lyapunov vectors. It should be 
emphasized that the wave-like structure characterizes the Lyapunov vector at any 
given instant and is not an average property. Our arguments depend essentially 
on the local and hyperbolic character of the interactions, in addition to translation 
invariance. The absence of translation invariance has been recognized to ruin the 
hydrodynamic modes [[7p. In translation invariant anharmonic chains, the absence 
of short time hyperbolicity seems to ruin the hydrodynamic modes [|8p. We present 
theoretical arguments for the existence of hydrodynamic modes in the simplified 
model, which are complemented by numerical verifications. The outline of the 
paper is as follows. In section ^ the hydrodynamic phenomenology is described 
in some more detail, the definition of the random matrix model is presented and 
motivated, and the results are stated. They are derived in sections |3|-|5]. Section || is 
devoted to numerical studies of the hydrodynamic properties of Lyapunov vectors. 
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2 Hydrodynamic behavior in tangent space 



The systems studied in [Q] consist (among others) of a large set of disks moving 
in a two dimensional box with periodic boundary conditions (torus geometry), 
with elastic scattering. In this case, the phase space is 4iV-dimensional where iV 
is the number of disks. The Lyapunov vectors have AN components which we 
label as 

(5x n , 5y n , 5p XjTl , 5p ytn ) 1 < n < N , 

with evident notation. To give them a geometrical meaning the components of the 
Lyapunov vectors are drawn in [Q] at the instantaneous position of the particles 
which carry a given specific index. That is, one constructs a vector field v(t, x) 
with values in R 4 , which are defined only at the instantaneous positions x n (t) of 
the particles, for example 

v x (t,x n (t)) = 8x n (t) , 

and similarly for the other components. 

The vector fields of the slow Lyapunov vectors, defined as the Lyapunov vec- 
tors with small corresponding Lyapunov exponents, are very well approximated 
by the long wavelength eigenmodes of a 'reverse wave equation' in the domain fh 

d^v(t,x) = -^V 2 v(t,x) (1) 

(note the unusual sign in front of V 2 ). That is, the vectors look like long wave- 
length waves with, say, n nodes in the x direction and m nodes in the y direction, 
and the corresponding Lyapunov exponent is proportional to 



1 /mr / n 

n\\lJ + {t 1 

Note that the translation modes — constant Lyapunov vectors with zero exponents, 
which are trivially present in any system with translation invariance, correspond 
to the special case m = n = 0. This phenomenology was observed in simulations 
with widely varying parameters, such as aspect ratio, density, and the shape of the 
particles 0,0]. 

The tangent flow of the molecular dynamics system can be written as 

*(*)-<*>(*)■ (2 > 

where the components 5x, 5p are column vectors with iV entries each of which 
is a vector in R 2 . The quantity G(t) is the action on the tangent space induced 
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by the flow $(£) of the dynamical system: If ip is the instantaneous state of 
the system then G(t) is given by G(t)f = D^(t)^ f, where $(t)(ip + ef) = 
&(t)ipo + eDQ{t)^ f + 0(e 2 ). The evolution operator of eq. (J2|) may be written 
formally as 



We proceed to construct a simplified model of the tangent dynamics, by making a 
series of modifications and assumptions about the nature of G and U. 

We first replace the hard-core interaction with a short range 'soft' potential. In 
that case G will have a block structure of the form 



where the symmetric N x N matrix A depends on the instantaneous positions of 
the particles and couples only nearest neighbors. Since the interactions are purely 
repulsive, the flow is hyperbolic, which implies that A should be taken positive. 

At this stage one may note that if the matrix A(t) in eq. <@> were replaced 
by the negative of the discrete Laplacian, the Lyapunov spectrum of G would 
be precisely that described in eq. ([[]). However, the matrices A(t) are in fact 
generated by chaotic dynamics, and therefore fluctuate rapidly. Furthermore, the 
particles in the gas rearrange in time, so that the positions of the non-zero elements 
in the matrix also evolve. In our study we concentrate on the first feature. That 
is, we show how tangent dynamics of the type (|]) result in slow hydrodynamic 
modes in spite of the fluctuations in A(t); the effects of particle rearrangement 
may in principle dealt with similarly, but need to be studied further. 

The above discussion allows us to conclude that the matrices A(t) should have 
non-zero elements only at those positions which are nonzero in the discrete Lapla- 
cian. Furthermore, momentum conservation implies that the sum of elements in 
any row and column of A must vanish. This specifies completely the matrix struc- 
ture of A, and it remains to model the time dependence of the off-diagonal non- 
zero elements of A. For this we invoke the hypothesis of strong chaos [JIJ, [Sp: The 
elements of A may be treated as independent random processes, with a correlation 
time r which is short with respect to other time scales of the system. It is com- 
monly found that this approximation yields results which are in good qualitative 
agreement with those of the actual tangent flow. 

With this in mind we model the evolution operator U by a product of indepen- 
dent random matrices S n 




(3) 




(4) 



u = Y[s, 



(5) 



n 



where 




(n-l)r 



(6) 
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During the time interval of length r, A, and therefore G, may be considered con- 
stant, so the simplest model for S would be S — ( ^ ^ \ However it is more 

convenient to correct this form by a second order term in r in order to preserve 
the symplectic property which holds for U. We thus arrive at our model [Q, [3j]: 

S= {rA lAt)' (7) 

The matrices A are independent, and their off-diagonal elements are independent 
and identically distributed. The actual probability distribution of the off-diagonal 
elements can be chosen arbitrarily, subject to the constraint of uniform hyperbol- 
icity, namely that the support of the distribution is strictly negative, and bounded 
away from zero. 

The model as defined above makes sense in any space dimension, but for the 
sake of simplicity we study it in one dimension. There it bears similarity to the 
tangent dynamics of an anharmonic chain. However, in the latter case the matrices 
S would be elliptic rather than hyperbolic. As we show below this is an essential 
ingredient in the mechanism for hydrodynamic modes, which are not present in 
the Lyapunov spectra of anharmonic chains [^]] . Unfortunately, we have not been 
able to find a model dynamics in real space whose tangent space dynamics would 
resemble that generated by the matrices of type ([7[). On the other hand our results 
do not use explicitly the dimensionality of the system and seem to be generalizable 
to higher dimensions. 

Since the individual matrices are symplectic, the Lyapunov exponents of (|]) 
come in pairs of equal absolute value and opposite signs. Translation and Galilei 
invariance imply the existence of two vanishing exponents. We concentrate our 
attention on the Lyapunov exponents Ajv_i and \n-2 of smallest positive value, 
and the corresponding Lyapunov vectors vn-i and v n-2- 

Before we go on, it is necessary to make precise what we mean by Lyapunov 
vectors. As is well-known, essentially all numerical methods for calculating the 
tangent space dynamics of the the Lyapunov vectors v n are defined as follows: 
One starts with an orthogonal matrix Q, mutiplies it from the left with the tangent 
matrix (in the present case S) and decomposes the result as SQ = Q'R, where 
Q' is orthogonal and R is upper triangular. This procedur is iterated to yield a 
sequence of Q t . The columns of the orthogonal matrices Q t are what we will call 
the Lyapunov vectors. The reader should note that these vectors are not the ones 
whose existence is proved in the multiplicative ergodic theorem. 

We can now state the main result of this paper. Existence of hydrodynamic 
Lyapunov modes: As the size N of the matrices tends to infinity the exponents 
Ajv-i and \n-2 as well as the vectors ujv-i andv N _ 2 are asymptotic to the expo- 
nents and vectors that would be obtained if the matrices A where replaced every - 
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where by the negative of the discrete Laplacian (properly rescaled). The statement 
holds for the Lyapunov vectors, which are random objects, in probability. 

The statement is spelled out only for two Lyapunov modes, which have nearly 
equal exponents, and where the deviation from hydrodynamic behavior is the 
smallest. However, as will become apparent from the arguments below, the result 
can be extended to a number of Lyapunov modes near the middle of the spectrum 
which is proportional to y/~N. Our numerical studies also indicate that this is in 
fact true. 

As already explained, the basic reason for the existence of the hydrodynamic 
modes is translation invariance. However, this general observation is not suf- 
ficient, and the actual proof is not trivial. It depends on the random nature of 
successive matrices, i.e., on strong chaos. Our strategy will be to show the ex- 
istence of hydrodynamic modes first in the spectrum of a single matrix of type 
A, i.e., a negative random Laplacian (sec. 0); then this will be used to show that 
such modes exist in the Lyapunov spectrum of non-symplectic products of type 
11(1 + A n ) in sec. |], which in turn will be used to show the same property for 
symplectic products in sec. 0. 



3 Spectral properties of a single matrix 



The matrix A defined in section [2] takes in one dimension the explicit form 

/ a\ + a 2 — a 2 • ■ ■ —a x \ 

-a 2 a 2 + a 3 -a 3 
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y — ai (J —dN a N + o,i ) 

where the a„ are positive identically distributed independent random numbers. 
The distribution of the a n is arbitrary, subject to the condition < a min < a < 
Omax with a min < a max , and normalized such that (a -1 ) = 1, for later convenience. 
The (•) always denote expectation with respect to the probability distribution of 
the a. We are not going to assume that the width of the distribution is small. 
The matrix A may be written as a product 

A = -dAd (9) 

where d and d are the discrete derivatives whose action on a vector v G Rw is 



















(8) 



(dv) n = V n - V n -i 



(dv) n = V n+1 - V n , 



(10) 



and A is a diagonal matrix with diagonal elements a n . (The indices are extended 
periodically so that a^+i = cti)- If all the a n were equal to one, — A would reduce 
to the discrete Laplacian matrix d 2 = dd. 
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We define the Fourier transform matrix F with elements 



? fcn = -L e ^*» s (id 



N 

which is a unitary transformation taking to C N , the subset of Cat (with stan- 
dard basis vectors e n ) consisting of vectors v for which V-k = v%, which is an N 
dimensional vector space over R. The components of A in the new basis are 

A H = (FAF^) k i = n*((a) 4,z + 5 fc _,)w , (12) 

where /i& = 1 — exp j^kj, and a is related to the Fourier transform of a 
considered as a vector in by 

a = N~ 1/2 F(a- (a)) . (13) 

The random variables dk are centered, and as sums of independent random num- 
bers their 'single-point' distribution is nearly Gaussian with variance 

~|2\ _ - («)" 



~a\') = , d4) 

so that they are typically small, of order 0{N^ 1 / 2 ). The joint distribution is not 
Gaussian. 

Note that // = 0, so that row and column of A are zero, with the translation 
vector e being trivially a zero eigenvector. We define the slow subspace 

^ = Span({eo,e 1 ,e_ 1 })nC 7V , (15) 

and its orthogonal complement, the fast subspace Vf. We will consider often below 
the block decomposition of A and other matrices into the fast and slow subspaces, 
e.g., 

a=( a : 1 • (i6) 



. A s{ A*. . 

Note that Vf contains slow as well as fast modes. 

The block A ss has small norm of order 0(N~ 2 ), and the off-diagonal blocks 
have norm of order 0{N^ 1 ). However, there are more specific properties of A 
which are needed to establish the existence of hydrodynamic eigenmodes. Con- 
sider the eigenvalue problem Av = \v. Letting v = du, and using the representa- 
tion @ gives an equation for u 

- d 2 u = \A~ x u . (17) 
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It is convenient to proceed by writing eq. ( |T7| ) in Fourier component form 

(|/i fc | 2 - X)U k = A h-qU q , (18) 

with the b k bearing a relation to a" 1 analogous to that between a k and a n , namely 

1 



kn{ 1, 



(19) 



Since < l/a^ n is 0(1) we find that the b k are 0{N~ 1 / 2 ) for the same reason 
that the a fc are. 

We claim that given a fixed m, and for N —>■ oo the system ( |?7| ) /las £wo 
linearly independent solutions u^ ±m \ X± m such that 



\u 



(±m) 



U 



(±m)i 



0(iV 



-1/2N 



A 



- 1 



I Mr, 



0(iV 



-l/2\ 



(20) 



(21) 



We justify the claim by showing that eqs. ( |20"1 ) and ( [21] ) are consistent with the 
eigenvalue equation (|TH|). For this we rewrite (|TH]) as 



~(m) 



A, 



iMfcl 2 - A r 



r,( m ) 

9 



(22) 



We assume that (PQ-pT|) hold; this implies that the sum over q in ( |22| ) is dominated 
by the two terms with q = ±m, that is, 



iMfcl _ I Mr, 



:(&*-„# + 6 fc+m «L m 2) , for ^ m. (23) 



On substituting this expression in the left-hand-side of ( |20| ) the sum over k is 
observed to be local, in the sense that it is dominated by terms with \k\ ~ m, 

where \fi k \ 2 ~ i^jf) ■ Since b k is 0(N~ 1 / 2 ), assumption ( |2"0"| ) is verified. On the 

gives 



other hand, using ( |20| ) in ( |2~2| ) for = m gives 

A 



I Mr, 



A, 



(24) 



Since 6 and 6 2m are 0(iV _1/2 ) it follows that A m /(|/i m | 2 - A m ) = 0(N 1/2 ) ver- 
ifying (|2~1|), which shows that (]2ll|-|2~T|) are indeed consistent with (JTSJ). 
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In terms of the original variables v, eq. ( |23| ) reads 



~(m) _ fJ'k^r 



V 



ih-mV^ + h+mV-m) , for \ k \ ^ m > ( 25 ) 



so that the norm of v±, the component of orthogonal to e ±m is small, 

1 _ , k 2 m 2 „ / 1 



(m) II 2 



In words, these eigenvectors are almost pure Fourier modes, i.e., eigenvectors of 
the discrete Laplacian. 

For further developments we also need to show these modes are the only ones 
with eigenvalues of order 0(N~ 2 ). 

This is established easily by noting that the sharp cutoff on the probability 
distribution of the a implies that every realization A satisfies the bounds 

- a m in<9 2 < A < -a max <9 2 , (26) 

and then by the minimax principle it follows that the pth eigenvalue of A is larger 
than a min times the pth eigenvalue of — d 2 (sorting the eigenvalues of both matrices 
in increasing order). 

The results of this section can be summarized using the decomposition of 
into slow and fast subspaces defined above. We have shown there exist small 
numbers e and A, and a number < a < 1, such that the matrix A can be block 
diagonalized, 

A = RDR T , RR T = 1 , D= ( D q d ) > ^ 

with the off-diagonal blocks bounded by ||-R s f||, ||-Rfs|| < e , an d me diagonal 
blocks obeying 

D { > A > a\ > D s > , (28) 

and furthermore 

||4»||<aA. (29) 

The orders of magnitude for A and e are e = 0(N~ 1 ^ 2 ) and A = 0(N~ 2 ), 
whereas a ~ 1/4. However, to keep the discussion reasonably general we are not 
going to use these specific values in our arguments below. Rather, we will make 
statements regarding arbitrary matrices which satisfy the conditions (|27|-|29|). 

Although this will not be used below, it is relevant to note that if we let m 
vary, the small parameter in ( |20[ ) and ( |2T| ) becomes m/N 1 / 2 . This means that we 
can expect a number of hydrodynamic eigenmodes which is proportional to \/~N ' . 
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Another way to see this is related to the study of the vibrations of one-dimensional 
disordered lattices which are modeled precisely by the eigenmodes of matrices of 
type (j^). There it is known that the localization length £ is proportional to A" 1 . 
Since A m ~ (2tt/N) 2 the localization length will reach N when m = 0(N 1 ^ 2 ). 
Thus, again, we only expect wave-like modes when m < 0(N 1 / 2 ). 



4 Products of matrices of the form 1 + r A 

In this section we use the properties derived in section |] to derive the existence 
of hydrodynamic modes in the Lyapunov spectrum of the product ]J n (l + rA n ), 
where the matrices A n are independent realizations of the random matrix defined 
in eq. (|]). Beside providing a step towards proving the existence of hydrodynamic 
modes in symplectic products, such a product may be regarded as the discrete 
approximation to a continuous tangent flow given by 

U = Texp J A(t)dt . (30) 

[Compare eqs. (Bh and (Q).] Although this does not correspond to the tangent flow 
of a mechanical system, it is nonetheless the simplest example where hydrody- 
namic Lyapunov modes can be expected. For convenience of further analysis we 
absorb r into the definition of A and change to Fourier basis once and for all, so 
the problem becomes that of a product 

n( i +^.)- (3d 

n 

Since the Lyapunov exponents of the slow part are expected to be smaller than 
the rest we aim at showing that the first N— 3 Lyapunov vectors span a subspace L f 
(of C at) which is almost orthogonal to V s in the sense that for any two unit vectors 
Uf E Lf and v s E V s one has |iif -v a \ 1. Although the subspace Lf changes after 
each step, we show below that the 'almost orthogonality' is propagated from step 
to step. 

To show this, we propose the following scheme. Take an arbitrary vector 
u E L f , whose components in Vf and are n f and u s respectively, normalized so 
that \\u{ || = I, and assume that u s is small. The action of 1 + A generates a new 
normalized vector u' by 

" = ir+iHii ■ <32) 

where [-] f is the projection onto the f -component We would like to show that u s 
remains small after repeated iteration of this process. 
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The block diagonalization ( |27| ) shows that the subspaces Vf and V s are indeed 
almost invariant under the transformation A. However, in trying to apply this 
fact to the Lyapunov vectors of the product ( |3~T| ) we immediately encounter the 
danger that the small perturbations may accumulate. The basic problem is that a 
vector in V B is contracted with respect to the 'slowest' direction in Vf by a factor 
of only 1 — 0(A) (as can be seen from the bounds on Df and D s ), whereas the 
perturbations which tilt a vector in Lf with respect to Vf are of order e, which is 
the typical size of the off-diagonal blocks [cf. eqs. (]27|)-(P9[)] and since we are 
interested in the case A <C e this contraction is not strong enough to overcome the 
perturbation. 

This order of magnitude argument can be made explicit by constructing a se- 
ries of matrices with the properties given by eqs. d27[)-(]29[), which take a vector 
in Vf and rotate it such that the outcome is a vector which has an angle with Vf of 
order 1. This counter-example is given in appendix [A|. 

An essential ingredient in the construction of this counter-example is that the 
matrix A has to be chosen specifically given u which in turn depends on for- 
mer realizations, in violation of the independence assumption. In other words, 
although such a 'bad' sequence is possible one naturally expects that this is an 
event with very low probability. Typically the perturbations to u s generated by 
the off-diagonal part of the matrices R do not have the same direction, and should 
serve to cancel one another. Therefore the statement one can hope to show is that 
in the sequence generated by iteration ofeq. (ffZ|), the probability that \\u s \\ > Ce 
for some fixed C is very small, as was shown for a similar example in [Q]. Here we 



will only prove the weaker statement that the variance 



2 ) is 0(e 2 ), and take 



that as an indication that the probabilistic statement is correct, since the behavior 
of higher moments can be treated in an analogous manner. 



To prove this statement we look at the s-component of eq. (32), 



ii a 



A si U{ + (1 + A ss )u s 
\(1 + A H )u{ + A is u s 



(33) 



The quantity (\\u' B \\ ) is a sum of three terms E\ + E 2 + E 3 : 



Ex 



|A,fWf|| 2 \ _ n / Agfut ■ (1 + A ss )u s \ _, / + A ss )u s \\ 2 ' 



En 



(34) 

where £ — ||(1 + A {i )u { + A {s u s \\. 

To bound these terms we first need a lower bound on the denominator £. Let 

Vf = Rf f Uf + Rf B u B , and define d by 



\\D iVi \\ =d\\v t \\ . 



(35) 
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Note that d can vary widely between 0(1) values and 0(A). But, using the lower 
bound on D { of (j28[), we see that 



||(1 + AHI 2 = INI 2 + 2v t ■ DfVf + \\D { Vi\\ 2 > (1 + 2A + d 2 )\\v { \\ 2 . (36) 
Expanding £ as 

£ = \\R&(1 + Df)v t + Rf s D s v s \\ , (37) 

and using the estimates ||i2f s D s || = 0(eX) and ||1 — R {i \\ = 0(e 2 ) (cf. eq. 
we get from (^6]) the desired lower bound on £: 



£ 2 >(l + 2X + d 2 )\\v i \\ 2 (l-0(s)). (38) 
We can now bound E\, E 2 and E 3 . First, we have 

IIAfMfH = H-Rsf-DfUfll + O(Xe) < sd\\v { \\ + O(Xe) . (39) 
Thus, neglecting higher order corrections in e, we get 

e 2 d 2 

Ei < . (40) 

l + 2A + rf 2 

The bound on the term E 2 makes essential use of the translation invariance. 
For this, we note that 

(41) 

||(l + A ff )n f || 2 ^ ^ 

transforms as a vector, that is, its kth component is multiplied by exp(z 2 ^/cx) un- 
der a relabeling of the coordinates n — > n + x. Therefore, because of translation 
invariance, the expectation value of ( pdj ) must remain invariant under such trans- 
formation, which means it must vanish. Since the denominator in E 2 is £ 2 (which 
also depends on u s ) and not || (1 + Aff)u{ || 2 , we need some gymnastics to exhibit 
the vanishing term. In order to see this we write 



Ei 



A s{ Uf ■ u s 



l + A ff W f | 



+ 



AfU{ ■ A ss u s 
{l + A H )u { \ 



A si u { ■ (1 + A ss )u s [2(1 + A {i )u { ■ A fs u s + \\A fs u s 
||(l + A ff K|| 2 £ 2 



(42) 



The first term in ( |42| ) vanishes because of translation invariance, as explained 
before. The second term is bounded by 



2a\e (\\u* 



(43) 
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and the dominant part of the third is 



(A S{U { ■ u s ){(l + A{ { ) Ui ■ Af s u s ) Ade 2 (\\u s \\ 2 ) 

+ A H ) Ul \\ 2 £ 2 l + 2A + d 2 ' ( ' 

The last term is bounded by 

l + 2a\ /„ 2 



g3< l + 2A + ^ (l|,/ - 



Collecting the bounds yields 



kll ) < 1 + 2 x + d 2 + 2aXe ■ (45) 



It appears from ( [431 ) mat although large perturbations are possible when d is 0(1), 
the contraction rate increases precisely enough to compensate this contribution. 
Thus if (||w s || 2 ) is 0(e 2 ) to start with, it will stay so indefinitely. 

In summary, assuming that the variance is indeed a measure of typical fluc- 
tuations, we have shown that, for N ^> 1, the subspace Lf spanned by the first 
N — 3 Lyapunov vectors of the product ( |3"T1 ) is, with very high probability, almost 
orthogonal to V^. This implies that the last three Lyapunov vectors (including the 
translation) remain approximately in V s . This means by definition that they are 
hydrodynamic, in the sense that they are well approximated by eigenvectors of 
the discrete Laplacian. Since the action of A on V s has two eigenvalues close to 
(2w/N) 2 as shown in the previous section, it follows as a corollary that the two 
smallest non-trivial Lyapunov exponents have approximately this value, so that 
they are also hydrodynamic. This completes the demonstration. 



5 Products of symplectic matrices 

We now turn to products of matrices of the form 

S =(rA 1+tm)' (46) 

We disregard the two translation modes in S for convenience and view the matri- 
ces S as (2N — 2) x (2iV — 2) matrices. Let us recall that since the matrices S 
are symplectic and hyperbolic, the Lyapunov exponents are non-zero and come in 
pairs of opposite signs. We concentrate on modes number N — 2 and N — 1 which 
are the smallest positive ones. 

We reduce the problem to an equivalent one to which the results of section |] 
can be applied directly. We denote by L + the subspace spanned by the first N — 1 
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Lyapunov vectors. It is spanned by a set of N — 1 independent vectors, which we 
display in the form of a (N — 1) x {2N — 2) matrix V. The N — 1 vectors can 
always be chosen in such a way that V is of the normal form 

V=(£), (47) 

where both blocks are (N — 1) x (JV — 1). Acting on V with 5 gives a spanning 
set of the image subspace L' + , 

SV=(_..\t (48) 



tA + (1 + t 2 A)V 
and changing basis to normal form gives V = (yi \ wnere 

V" = r A + — ^— . (49) 

1 + TK 

A convenient property of this matrix dynamical system is that if V is symmetric 
to begin with, it stays so as a consequence of the symplectic property of S [|10l]. 

By definition, any vector v E L + has a block representation v = (u, Vu). 
From eq. ( ^S[ ) it follows that its image is ( u', V'vl ), where 



it' 



(1 + tV)u . (50) 



Hence, the first N — 1 Lyapunov modes of the products of the S are the same as 
those of the product Yl n (l + rV n ) where the matrices V n are evolving according 
to eq. ®: V n+1 = rA n + V n /(1 + rV n ). 

In view of this equivalence, it suffices to show that the matrix V has the prop- 
erties formulated in eqs. (p7|-p9|) and to apply the results of section |]. First note 
if A n = —d 2 (minus the discrete Laplacian) for all n, then the V n converge to 
f(—d 2 ), where 




tx x 2 



f{x) = — + Jx + 



is the larger root of the quadratic equation f(x) = tx + t+^jm ■ For small 

x > this is close to x 1 ^ 2 , and therefore we assume that V has a representa- 
tion of the type given by eqs. (|2"7|-|2"9D, with e = 0{N~ 1 / 2 ) as before and A is now 



f(4:%/N 2 ) = OiN' 1 ). The aim is to show that this property is carried on to V. 

In order to avoid the necessity of presenting even more technical details, we 
present the argument for the case where the slow subspace contains a single mode, 
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rather than a pair of nearly degenerate modes. Since V is symmetric its smallest 
eigenvalue is given by 

V 

Ay/ = min u ■ (tA + - —)u . (51) 

IMI=i 1 + tV 

It follows from our assumptions that there exist (normalized) eigenvectors of V 
and A: 

Ae A = Ai(l + c A e)e A , Ve v = /(Ai)(l + c v e)e v , (52) 

where Ai is the smallest positive eigenvalue of — d 2 and e A and ey are close to e\, 
the corresponding eigenvector. The variational principle then gives immediately a 
lower bound on Ay/, 

Ay, >/(A a )(l + ec v ,), (53) 

for some cy/ between c A and cy. To get an upper bound on Ay/ recall that it was 
shown in sec. 0eq. ( |2~5l ) that the component of eA is of order ek' 1 , and note that 
the bound d26[ ) implies 



/(-<w9 2 ) < V < /(-a max «9 2 ) . (54) 
We now use u = e A in eq. ( |5T] ) and get 

Ay/ < Ai(l + c A e) + e A ■ f(-a max d 2 )e A < f(Xi)(l + ce) , (55) 



for some constant c. Eqs. ( p3[ ) and ( p5D establish the desired property of the eigen- 
values. 

The corresponding eigenvector ey is the one which minimizes eq. (|5T|). Be- 
cause of the minimax principle applied to A and V, letting u = e\ + w with 
w • e\ = and u» small, the quadratic form u ■ Au may be approximated by 

u ■ Au ~ (w — tUyi) • 74(w — Wi) + Aa , 

and similarly 

M • Vu ~ (w — %) • — lOy) + Ay . 

Therefore, in order to find w V i we need to minimize the quadratic form 

V 

t(w - w A ) ■ A(w - w A ) + (w- w v ) ■ — — —(w - w v ) . 

1 + TV 

The minimum occurs at 

wv> = (1 + B)-\Bw A + w v ) , (56) 

where 

1 + TV a 

B = r^^A. 
15 



We can use again the bounds ( [26] ) and (Q) to show that 



b m \n9{-d 2 ) < B < b max g(-d 2 ) , 

for some positive numbers b min , b max and a positive function g, and thus bound 
the components of wv>, 

u \ | (wv)k + b nmx g(\nk\ 2 )(w A )k 
\ w v')k < TTT 7\ — m ' 

1 + 0min#(|/-tA:r) 

where \nk\ 2 is the A;th eigenvalue of — d 2 (see sec ^]). This shows that (wy')k ~ 
ek~ l . This completes the demonstration of the desired properties of V, and, on 
applying the results of section ^ the existence of hydrodynamic Lyapunov modes 
in the symplectic case. 



6 Numerical tests 

The purpose of this section is to verify numerically some of the statements given 
above, and on to further study numerically the dependence of hydrodynamic be- 
havior of several Lyapunov modes as a function of noise level as well as system 
size. 

The simplest system we discuss is a product Y\ n A n of independent matrices 
of the form (|8|). Since the relative gap between the first two non-zero eigenvalues 
is 0(1) [see eqs. (|Z7|-|Z9|)], the contraction in this case is strong, and the potential 



problems of the accumulation of errors as discussed in sec. and appendix |A| 
are absent. Nevertheless, even in this case, there are some qualitative differences 
between the behavior of the Lyapunov modes, and the corresponding eigenmodes 
of a single matrix. 

We quantify the degree of hydrodynamic behavior in the Lyapunov modes as 
follows. For the Lyapunov vectors we computed the residuals r^, that is, the 
norm of the orthogonal complement 

n = Ik - (vi ■ e k )e k - (vi ■ e_ fc )e_ fc || 

where k — k(i) is the wave vector associated with vector i. (For example, k — 1 
for the vectors Ujv-i an d v jv-2 discussed above.) In fact, to get more precise results 
we subtracted from Vi all the components with lower-lying k: 

r i = Ik - (( v 'i ' e k) e k - ( v i ■ e_ fc )e_ fc )|| 

k<k(i) 

(The results are not very different for the two definitions of ?";.) 
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Deviation from Eigenspace, 1 — A 
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Figure 1: Data collapse of the residuals of the Lyapunov vectors of a product of 
random Laplacians. The curves show (rt/cr) 2 as a function of the wave number 
k/N for several values of iV and a. The data points are averages of 10 realizations, 
averaged also within each pair which corresponds to the same \k\. As expected, 
the range of collapse increases with the system size. The random variables are 
chosen as a = 0.5 ± a with probability 1/2. 
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Fig. [T| presents these residuals for systems with different sizes of the matrices 
and different values of the noise variance a. The vertical axis measures (r k /a) 2 , 
and the horizontal axis gives k/N. The approximate collapse of the graphs for 
small k implies that the dependence of the residuals on system size N and noise 
strength (variance) a) is given by the scaling form 

The behavior of /i for small x is approximately fi(x) = 0(y/x), which implies 
that for a fixed k 

a 




'N 

the same dependence as in the residuals of the eigenvectors of a single matrix. 
However, the dependence of the residuals as a function of k is r k ~ yk, slower 
the linear dependence on k in the case of a single matrix. 

For the Lyapunov exponents A» we measure the relative deviation 5k from the 
respective eigenvalues |/i fc | 2 of the discrete Laplacians (see sec. Q), 

exp(A fc ) 



2 



- 1 




The results for the deviations Si of the Lyapunov exponents are displayed in fig. |2|, 
where 5k/ a 2 is plotted against k/N. The data collapse implies that 

$k = 

The function /^(^O is approximately linear for small x which implies for k <C N: 

This is not unreasonable, since the Lyapunov exponents, unlike the eigenvalues of 
a single matrix, are given as a result of an averaging process. 

We next present a similar analysis for the product Y\ n (l + A n ) which was 
considered in section |]. The results for the residuals of the Lyapunov vectors and 
the deviations of the exponents are presented in figures |3] and |] respectively. The 
scaling form for the residuals is in this case 

r k ■ 

where f^x) ~ x for small x. Thus, in this case the residual for a fixed k decreases 
as 




a 

rk 



N ' 
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Difference of Liap. exponents, — A 
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Figure 2: Data collapse of the relative deviations of the Lyapunov exponents of a 
product of random Laplacians. The curves show 5 k /a 2 as a function of the wave 
number k for several values of N and a. In fact the data collapse is perfect for 
fixed a. 
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that is, faster than the N~ 1 / 2 decrease in the residuals of the eigenvectors of a 
single matrix. The analysis presented in section [| is too general to capture this 
behavior. 

The relative deviations of the exponents scale in this case as 

h = <y 2 h (j^j , 

and fi(x) ~ x 2 for small x, so that as in the product of random Laplacians, the 
relative size of r k and 5 k is 5 k ~ r\ . 



A Counter-example 

We want to show here that an 'unfortunate' choice of rotations can move the sys- 
tem out of the region where the Lyapunov vectors remain essentially aligned with 
the eigendirections of the Laplacian. The issue here is that, on one hand, the cones 
in which these vectors lie are slightly contracted and on the other hand slightly 
turned. The 'counter-example' shows that the turning wins over the contraction. 

Let Vf = Spanjei, e 2 } and V s = Span{e 3 }. Suppose that L f contains a 
vector with block representation u = (u f , u B ), with u s > 0, normalized such that 
||«f || = 1, and let t>f span the orthogonal complement to Uf in Vf . We construct the 



matrix A by giving the components in the representation (27). 



L> s = aX, and Dr is such that 



-DfWf = (A + £ 2 )u { + evf 



(58) 



The image of u is 



(l + A)( Ui ) = ( il+ ( ^ £ ^- 2 ^ )u jV Vi ) ■ (59) 
\u s J \ (1 + ocX - e )u s + e A ) 

After normalizing the / component to 1 , the s component becomes 



l + («-l)A + e 2 (3/2-M s ) u s + e 2 . (60) 



Evidently, even if u s = initially, by choosing A as above, u s can be increased to 
an 0(1) value (as e, A -> 0) if A = 0(e 2 ). 
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Deviation from Eigenspace, 1 — A 




Figure 3: Same as fig. |T] for the product H n (l + A n ), except that the vertical axis 
measures r k /a. 
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Difference of Liap. exponents, 1 — A 




Figure 4: Same as fig. ||for the product H n (l + A n ), except that the vertical axis 
measures yfd^/a. 



22 



Acknowledgments: We have profited from very useful discussions with S. 
Lepri, C. Liverani, Z. Olami, A. Politi, and L.-S. Young. But we are particularly 
indebted to H. Posch for having provided and discussed with us his beautiful nu- 
merical experiments. This work was partially supported by the Fonds National 
Suisse, and part of it was done in the pleasant atmosphere of the ESI in Vienna 
and at the Istituto Nazionale di Ottica in Florence. 

References 

[1] Paladin, G. and Vulpiani, A.: J. Phys. A19, 1881 (1986). 

[2] Livi, R., Politi, A., Ruffo, S., and Vulpiani, A.: J. Stat. Phys. 46, 197 (1987). 

[3] Eckmann J.-P. and Wayne, C.E.: J. Stat. Phys. 50, 853 (1988). 

[4] Dellago, Ch., Posch, H.A., and Hoover, W.G.: Phys. Rev. E53, 1485 (1996). 

[5] Posch, H.A. and Hirschl, R.: Simulation of billiards and of hard-body fluids, 
preprint (1999). 

[6] Forster, D.: Hydrodynamic Fluctuations, Broken Symmetry, and Correlation 
Fluctuations, Benjamin Cummings, London, (1983). 

[7] Lepri, S., Livi, R., and Politi, A.: Europhys. Lett. 43, 271 (1998). 

[8] Politi, A.: private communication. 

[9] Young, L.S.: Erg. Th. Dyn. Syst. 6, 627 (1986). 

[10] Partovi, M.H.: Phys. Rev. Lett. 82, 3424 (1999). 



23 



